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Abstract 

We report on quantum Monte Carlo calculations of the ground and low-lying excited states of 
A = 9, 10 nuclei using realistic Hamiltonians containing the Argonne v\% two-nucleon potential 
alone or with one of several three-nucleon potentials, including Urbana IX and three of the new 
Illinois models. The calculations begin with correlated many-body wave functions that have an 
a-like core and multiple p-shell nucleons, LS*-coupled to the appropriate (J W ;T) quantum numbers 
for the state of interest. After optimization, these variational trial functions are used as input to 
a Green's function Monte Carlo calculation of the energy, using a constrained path algorithm. We 
find that the Hamiltonians that include Illinois three-nucleon potentials reproduce ten states in 
9 Li, 9 Be, 10 Be, and 10 B with an rms deviation as little as 900 keV. In particular, we obtain the 
correct 3 + ground state for 10 B, whereas the Argonne v±s alone or with Urbana IX predicts a 1 + 
ground state. In addition, we calculate isovector and isotensor energy differences, electromagnetic 
moments, and one- and two-body density distributions. 

PACS numbers: 21.10.-k, 21.45. +v, 21.60.Ka, 27.20. +n 



I. INTRODUCTION 



In a series of papers |], ||, 0, we have reported quantum Monte Carlo (QMC) calculations 
of ground and low-lying excited state energies in A < 8 nuclei for realistic nuclear Hamiltoni- 
ans. These calculations employed the Argonne v is (AV18) two-nucleon potential 0] and the 
Urbana IX (UIX) three-nucleon potential |IJ, and are accurate to ~ 1 — 2% of the binding 
energy for light p-shell nuclei. More recently, we have used the quantum Monte Carlo cal- 
culations to help construct a series of new and improved pion-exchange-based three-nucleon 
potentials, designated the Illinois models [||]. The five Illinois models (IL1-IL5), when used 
in conjunction with AV18, each reproduce the experimental energies of 17 narrow states 
in A < 8 nuclei with an rms deviation of ~ 400 keV. This contrasts with a 2.3 MeV rms 
deviation for the AV18/UIX Hamiltonian, and 7.7 MeV for AV18 alone. 

In this paper, we report the extension of our QMC calculations to A = 9, 10 nuclei. The 
QMC methods include both variational (VMC) and Green's function Monte Carlo (GFMC) 
methods. The VMC method is used to construct a variational wave function as a product 
of two and three-body correlation operators acting on a independent-particle wave function 
that has an a-like core and multiple p-shell nucleons, L.S'-coupled to the appropriate ( J n ; T) 
quantum numbers for the state of interest. Monte Carlo evaluation of the energy expectation 
value is used to optimize the trial function, particularly the mix of independent-particle 
wave function components. The GFMC method starts from this trial function and makes 
a Euclidean propagation that converges to the lowest energy for a state of these quantum 
numbers. A constrained path algorithm is crucial for keeping the fermion sign problem 
under control. At present the GFMC method is used to calculate only the lowest state of 
given (J n ;T). 

We have calculated ten stable or very narrow natural-parity states in the A = 9, 10 
nuclei 9 Li, 9 Be, 10 Be, and 10 B that are experimentally well known. We use three of the 
new Illinois models (IL2, IL3, IL4) in conjunction with AV18 and obtain rms deviations 
from the experimental energies of the states of 900, 1100, and 1700 keV, respectively. We 
have also calculated most of these states with the AV18 and AV18/UIX Hamiltonians for 
comparison. The most intriguing result we find is that the new AV18/IL2-IL4 models all 
correctly predict a 3 + ground state for 10 B, but the older models wrongly predict a 1 + ground 
state. In addition, we have made GFMC calculations of six other states that are expected 
within the p-shell formulation of these nuclei and of the 9 He and 10 He ground states; these 
states either have much larger widths or are not clearly identified by experiment. 

We review briefly the experimental status of the ground and low-lying excited states in 
the A = 9, 10 nuclei in Sec. II. The Hamiltonians are described in Sec. Ill and the QMC 
calculations in Sec. IV. Most of this material has been discussed in detail in Refs. 0, 
The only major new technical development is the automation of the construction of the 
independent-particle portion of the variational trial wave functions that serve as starting 
points for the GFMC calculations. Energy results are given in Sec. V, while electromagnetic 
moments and density distributions are shown in Sec. VI. We present our conclusions in 
Sec. VII. 

II. EXPERIMENTAL STATUS 

The experimental status of A = 9, 10 nuclei is illustrated in Figs. [I] and |2|, where we 
show the ground states and most low-lying natural-parity states whose spin assignments are 
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reasonably certain || |7| . We also show some additional narrow states in 9 Li and 10 Be whose 
spins have not been determined experimentally; reasonable guesses are given in parentheses, 
based on standard shell model studies and the present calculations. 

The ground state of 9 Be is an absolutely stable (J 77 ; T) = (| ; |) state. The first exci- 
tation is a (| + ; |) state (not shown in Fig. |l|) with a width of 217 keV which occurs at the 
threshold for breakup into 8 Be + n. This unnatural parity state is beyond the scope of the 
present paper; we will report calculations of such intruder states in the future. The second 
excitation is a narrow (| ; |) state at 2.429 MeV, with a width < 1 keV. Within the p-shell 

formulation, there can also be (| ; |), (| ; |), and (| ; |) states. Experimentally, the first 
two are observed at 2.78 and 6.38 MeV excitation but both are quite broad (~ 1 MeV), 
while no state with (| ; |) character has been identified. We evaluate all these (J n ;T) cases 
in GFMC, treating the states as if they were particle stable; this should be adequate for 
narrow states, but may be less satisfactory for broad states. (Several additional unnatural 
parity states and second excited states of given ( J n ; T) are observed above 3 MeV, but are 
not shown in the figure.) The matrix elements of the electromagnetic and strong charge- 
independence-breaking terms in the Hamiltonian are evaluated perturbatively to infer the 
energies of the narrow T — ~ isobaric analog states in 9 B. 

The ground state of 9 Li is a particle stable (| ; |) state that decays by (3~ emission to 
9 Be with a half-life of 178 ms. The first excited state of 9 Li at 2.691 MeV is believed to 
be a (| ; |) state; it is below the threshold for breakup into 8 Li + n and decays only by 
7 emission. We calculate both these states in 9 Li, and again perturbatively evaluate the 
energies of the isobaric analogs in 9 Be, 9 B, and 9 C. Two additional narrow states above the 
breakup threshold have been observed in 9 Li, but without firm spin-parity identification: a 
60 keV wide state at 4.296 MeV and a 40 keV wide state at 6.43 MeV. Within the p-shell, 
these should be (| ; |) and (| ; |) states; our GFMC evaluations of such states line up very 
well with the experimental observations, suggesting these spin assignments may indeed be 
correct. 

The ground state of 10 B is an absolutely stable (3 + ;0) state. The first threshold for 
breakup of 10 B is into 6 Li+<x Between the ground state and breakup threshold are two 
(l+;0) excited states at 0.718 MeV and 2.154 MeV, one (2+;0) state at 3.587 MeV, and the 
(0 + ;l) isobaric analog of 10 Be at 1.740 MeV. Many additional states are known above the 
6 Li+a threshold; in Fig. |2| we show only the (2 + ;l) isobaric analog at 5.164 MeV and the 
(4+;0) state at 6.025 MeV. We calculate the (3 + ;0) ground state and first (l + ;0), (2+;0), and 
(4 + ;0) excitations by GFMC. In the p-shell there can also be (0 + ;0) and (5 + ;0) states, but 
the former is of low spatial symmetry while the latter has high angular momentum, leading 
us to expect both of them to be quite high in excitation energy; no experimental observation 
of either has been made. 

The nucleus 10 Be is particle-stable but decays by f3~ emission with the very long half-life 
of 1.51 x 10 6 years. As is typical for an even-even nucleus, the ground state is a (0 + ;l) state, 
and there is a well-separated (2 + ;l) excited state at 3.368 MeV. These are the two primary 
10 Be states we evaluate in GFMC. In addition, the tabulation lists a state at 9.4 MeV as 
a possible (2 + ;l) state, but a more recent (t, 3 He) experiment makes a (3 + ;l) assignment 
more likely || |9j. We have made one computation for the (3 + ;l) level using AV18/IL2 and 
get reasonable agreement with this energy. Other particle-stable states below the threshold 
for breakup into 9 Be+n include second (0 + ;l) and (2 + ;l) states, also shown in Fig. ||, and 
the first particle-stable (1~;1) and (2~;1) intruder states (not shown). The second excited 
states are not evaluated at present in GFMC, while the intruder states will be reported on 
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in future work. Again, the isobaric analog states in 10 B and 10 C are evaluated perturbatively 
from 10 Be. 

In addition to these nuclei, we have also made calculations of the expected lowest natural- 
parity states in 9 He and 10 He. These are a (~ ; ~) resonance which is observed at ~ 1.2 MeV 
above the threshold for breakup into 8 He+n, and a (0 + ; 3) resonance at « 1.1 MeV above 
the 8 He+n+n threshold |7|, [nj. These resonances are observed to be reasonably narrow, 
with widths of 100 and < 300 keV, respectively. There are recent experimental reports of 
a (| + ; |) resonance near threshold in 9 He [O, which we have not attempted to calculate; 
such a state is likely to be very broad and really needs to be treated as a scattering state. 
We also have not attempted to calculate the resonant ground state of 10 Li, which is broad 
and still has some experimental uncertainty. 

III. HAMILTON! AN 

The Hamiltonian includes nonrelativistic one-body kinetic energy, the AV18 two-nucleon 
potential M and either the UIX |TJ or one of the Illinois 0] three-nucleon potentials: 

H = Y f K i + y £v ij + ]T V ijk . (1) 

i i<j i<j<k 

The kinetic energy operator is predominantly charge-independent (CI), but has a small 
charge-symmetry breaking (CSB) component due to the difference between proton and neu- 
tron masses. The AV18 is one of a class of highly accurate NN potentials that fit both pp 
and np scattering data up to 350 MeV with a x 2 /datum « 1. It can be written as a sum of 
electromagnetic and one-pion-exchange terms and a shorter-range phenomeno logical part: 

%=4 + v S + <- ( 2 ) 

The electromagnetic terms include one- and two-photon-exchange Coulomb interaction, vac- 
uum polarization, Darwin-Foldy, and magnetic moment terms, all with appropriate proton 
and neutron form factors. The one-pion-exchange part of the potential includes the small 
charge- dependent (CD) terms due to the difference in neutral and charged pion masses. 
The shorter-range part has about 40 parameters which are adjusted to fit the pp and np 
scattering data, the deuteron binding energy, and also the nn scattering length. 

The one-pion-exchange and the remaining phenomenological part of the potential can be 
written as a sum of eighteen operators, 

«5 + v 5 = £ «pfa)a& • (3) 

p=l,18 

The first fourteen CI operators are 

Of= 1M = [1, (Ji ■ aj, Sij, L ■ S, L 2 , LVi ■ <Tj, (L ■ S) 2 ] ® [1, r t ■ Tj ] , (4) 
while the last four, 

are the strong interaction CD and CSB terms. 
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The three-nucleon potentials from the Urbana series of models |T2| contain a long-range, 



two-pion-exchange, P-wave term and a short-range phenomenological piece: 

vSu = + ■ (6) 

The UIX model has the strengths of these two terms adjusted to reproduce the binding 
energy of 3 H, in GFMC calculations, and to give a reasonable saturation density in nuclear 



matter, in variational chain summation calculations [jT3|| , when used with AV18. The Illinois 



models add a two-pion-exchange S-wave term and a three-pion-ring term: 

VH = V$ p + V% s + V^ R + Vi % . (7) 

The two-pion-exchange S-wave term is required by chiral symmetry, but in practice its small 
energy contribution makes it hard to distinguish from the dominant P-wave term. However, 
the three-pion-ring term, while it is smaller than the two-pion-exchange P-wave term, has 
a distinctly different isospin dependence, which is crucial for being able to fit the variety of 
light p-shell energy levels studied in Ref. ||. In the Illinois models, the operator structure 
and radial forms were taken from standard meson-exchange theory, but the overall strengths 
of the four terms, and one cutoff factor in the radial dependence, were adjusted to obtain 
best fits to the energies of 17 narrow states in 3 < A < 8 nuclei. In practice, at most three 
parameters at a time could be uniquely determined from the energy calculations, so five 
different models were constructed in which different subsets of the parameters were fixed by 
external considerations, while the remaining ones were adjusted. 

The CD and CSB terms in H are fairly weak, so we can treat them conveniently as a 
first-order perturbation and use a wave function of good isospin, which is significantly more 
compact. Also, direct GFMC calculations with the spin-dependent terms that involve the 
square of the orbital angular momentum operator can have large statistical fluctuations, as 
discussed in Ref. 0. Thus it is useful to define || a simpler isoscalar interaction, AV8', which 
contains only the eight operators [1, Cf • (Tj, SV,-, L • S] ® [1, Tj • Tj] and an isoscalar Coulomb 
interaction. These eight operators are chosen such that AV8' reproduces the CI part of the 
full AV18 interaction in all S and P waves as well as the deuteron. The AV8' interaction 
(without the Coulomb term) was recently used in a benchmark test of seven different many- 
body methods for solving the four-nucleon bound state, with excellent agreement between 
GFMC and the other calculations 1131 . 



IV. QUANTUM MONTE CARLO 
A. Variational Monte Carlo 

We first construct a variational wave function for the state of interest and then optimize 
it by minimizing the energy expectation value as computed by Metropolis Monte Carlo 
integration. The variational wave function for the nuclei studied here has the form 

l*v>=[l+ E U^U^ + Ui^j). (8) 

i<j<k i<j 

The Uij, and U™ 1 are noncommuting two- and three-nucleon correlation operators, and 
S is a symmetrization operator. The includes spin, isospin, and tensor terms induced 
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by the two-nucleon potential, while the Uf^ 1 reflects the structure of the dominant parts 
of the three-nucleon interaction. This trial function has the advantage of being efficient to 
evaluate while including the bulk of the correlation effects. A more sophisticated variational 
function can be constructed by including two-body spin-orbit correlations and additional 
three-body correlations, as discussed in Ref. ||, but the time to compute these extra terms 
is significant, while the gain in the variational energy is relatively small. Studies have shown 
that the GFMC algorithm easily corrects for the omission of these terms 0. 

The two-body correlations are generated by the solution of coupled differential equations 
with embedded variational parameters []T5|. We have found that the parameters optimized 



for the a-particle are near optimal for use in the light p-shell nuclei. Likewise, the best 
parameters for the three-body correlations are remarkably constant for different s- and p- 
shell nuclei, so they have not been changed significantly from the previous A < 8 work 

For the p-shell nuclei studied here, the totally antisymmetric Jastrow wave function, 
starts with a sum over independent-particle terms, <& A , that have 4 nucleons in an a-like core 
and (A — 4) nucleons in p-shell orbitals. We use LS coupling to obtain the desired JM value 



of a given state, as suggested by standard shell-model studies fll6fl . We also need to specify 
the spatial symmetry [n] of the angular momentum coupling of the p-shell nucleons |17| . 
Different possible LS[n] combinations lead to multiple components in the Jastrow wave 
function. This independent-particle basis is acted on by products of central pair and triplet 
correlation functions, which depend upon the shells (s or p) occupied by the particles and 
on the LS[n] coupling: 

= a { n n fsMi) n /*(*•«) 

i<j<k<4 i<j<4 k<4<l<A 

E (PLS[n] II 4pW) \®A(LSin]JMTT 3 ) 12U:5 ... A )) } . (9) 

LS[n] 4<Km<A 

The operator A indicates an antisymmetric sum over all possible partitions of the A particles 
into 4 s-shell and (A — 4) p-shell ones. The pair correlation for particles within the s-shell, 
f ss (r), is the optimal correlation for the a-particle. The f sp (r) is similar to the f ss (r) at 
short range, but with a long-range tail going to a constant « 1; this allows the wave function 
to develop a cluster structure like a + d in 6 Li or a + a in 8 Be at large cluster separations. 
The fj$(r) depends on the nucleus and particular independent-particle channel, e.g., in the 
case of 6 Li or 8 Be, it is similar to the optimal deuteron or alpha correlations. 
The LS[n] components of the independent-particle wave function are given by: 

\$ A (LS[n] JMTT 3 ) 12Ut5 .„ A ) = |$ Q (0000) 1234 II # SW (^rf) 

4<KA 

[ II Y lmi (n Ql )] LMli[n] ®[ II Xi(lrn s )] SMs \ (g>[ II ui(h 3 )]TT a ), (10) 

4<1<A 4<1<A Z ) JM 4<KA Z 

where 

$ a (0000) = ^(j9tp|^Tn|) (11) 

is the a-core independent-particle wave function. The (pp S (R a i) are p-wave solutions of a 
particle of reduced mass |m n in an effective a-N potential: 

V aN {r) = V%f( r ) + V° N (r) . (12) 
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The <pp S are functions of the distance between the center of mass of the a core (which 
contains particles 1-4 in this partition) and nucleon /, and again may be different for different 
LS[n] components. For each state considered in the present work, we have used bound-state 
asymptotic conditions for the 0p 5 , even if the state is particle unstable. The Woods-Saxon 
potential 

V^ S (r) = V p LS [l + eM r —^T\ (13) 

tip 

has variational parameters V p s , R p , and a p , while the Coulomb potential is obtained by 
folding over nuclear and nucleon form factors. The wave function is translationally invariant, 
hence there is no spurious center of mass motion. 

A major technical advance in the present work is the automatic generation of the 
independent-particle wave function $^ with the appropriate spatial symmetries discussed 
above. We describe here the construction of the spatial symmetry [n] . One can use different 
coupling schemes to form the spatial (or spin and isospin) part of the wave function with 
good quantum numbers. We use simple sequential coupling in which the spatial part of the 
wave function is written as 
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LM L n 



a7j 



h&7 



(14) 

LMr 



and a similar construction is used for the spin and isospin part. The functions having the 
same L and Ml but different intermediate quantum numbers, labeled by /i, are orthogonal 
and form a complete set of eigenstates of L 2 and L z . The permutation operators of the 
valence particles I — 5,6, A (k = 1, (A — 4)!), commute with L 2 and L z so that the 
above functions form a basis for the representation of the symmetric group as well: 

PkQhM Li i = J~] U\ix{Pk)^LM,.\ ■ (15) 
A 

The permutation symmetry is conveniently depicted by using "Young diagrams" , consist- 
ing of N adjoining square boxes with rows numbered numerically downward, and columns 
rightward; there may not be more rows in column % than in column i — nor columns in row 
j than in row j — 1. Each Young diagram corresponds to a representation of the permutation 
group. The basis functions defining a given representation can be labeled by using a Young 
tableau, which is an arrangement of the numbers 1, 2, TV in the Young diagram, such that 
numbers always increase along all rows and down all columns. 

We use the "Young operators" to construct a basis function that has the symmetry 
properties of a given Young tableau. The Young operator, Y, is a product of a symmetrizer 
R and an antisymmetrizer Q: Y = QR. The operator R symmetrizes all particle indices 
which are in the same row and Q antisymmetrizes all particles in the same column. Both 
operators are constructed as a combination of the permutation operators 

To construct the spatial part of the wave function belonging to a given representation of 
the permutation group, we first have to draw the Young diagram and insert the numbers 
1,2, ...,N into the pattern in any order to give a Young tableau. The number of different 
tableaux Ny is the dimension of the representation. Then we prepare the Young operators 
Yi, ...,Y Ny corresponding to the Young tableau and calculate the matrix elements of these 
operators with the basis functions 0lm l h- We then have 

YkQhM L p, = ^ ( ®LMt . ix 1 Y k 1 9 LM r . A ) @LM r . A ■ (16) 
A 



7 



Let us denote the successively coupled spin (isospin) functions by ujsm s uu (<^tt 3 j/)- To form 
an antisymmetric wave function for the iV particle system one has to multiply the basis 
functions of the space part of a given Young tableau, Y^9lm l h, by the basis functions of the 
spin-isospin part belonging to the conjugated [n] Young diagram (obtained by reversal of 
the roles of rows and columns) Yk[n]uJsM s u^TT 3 u' and sum over all possible tableaux. Thus, 
Eq.fllOD becomes: 



\<5> A (LS[n]JMTT 3 ) 12U:5 ... A ) = |$ a (0000) 1234 II $ S[n] (R*i) 

4<1<A 

^PkYk[n]Y k [n][6 L ^ Sv } JM u TT , iV ,) , (17) 



k 

where p k is the parity of the permutation of the numbers (starting from the top and going 
left to right in each row) in the Young diagram. The Young operator Y&[n] (Y/Jn]) acts on 
the spatial (spin-isospin) functions of the I = 5, 6,..., A valence particles. Any choice of \i 
and v, v' generates the same wave function. 

The different possible LS[n] contributions to A = 9, 10 nuclei are given in Table |. After 
other parameters in the trial function have been optimized, a series of energy evaluations 
are made in which the (3ls[u] of Eq.(§) are different in the left- and right-hand-side wave 
functions to obtain the diagonal and off-diagonal matrix elements of the Hamiltonian and the 
corresponding normalizations and overlaps. The resulting NxN matrices are diagonalized to 
find the Plsm eigenvectors, using generalized eigenvalue routines because the correlated \lV 
are not orthogonal. This allows us to project out not only the ground-state trial functions, 
but also excited-state trial functions of the same ( J n ; T) quantum numbers. In our present 
studies, we have carried out our A=9 diagonalization in a complete p-shell basis, but for 10 B 
and 10 Be we have limited ourselves to the three highest spatial symmetries, i.e., [42], [33], 
and [411]. The diagonalization is carried out for the AV18/UIX Hamiltonian; the (3ls[u] 
amplitudes should not be significantly different for the other models. Additional spatial 
symmetries involving particle excitations out of the p-shell are built up in the full trial 
function by means of the tensor correlations contained in the Uij and U™ 1 of Eq.(||). 

Thus in 9 Be, the Jastrow wave function for the (| ; |) ground state is constructed from 
thirteen amplitudes, and a 13 x 13 diagonalization is performed to find the optimal mixing. 
The Plsm values for this and other 9 Be states, including some second excited states, are 
given in Table [TT], while the values for various states in 9 Li are given in Table JTJ. In the case 
of 9 Be, the [41] symmetry states dominate; addition of the lower spatial symmetries improves 
the energies by typically 0.25 MeV. The additional states do not significantly alter the rms 
radii or electromagnetic moments. By contrast, the leading [32] symmetry in 9 Li is not so 
dominant; addition of lower spatial symmetries gives a significant 1-2 MeV improvement to 
some of the energies. Electromagnetic moments can also be significantly shifted. In general, 
the dominant A = 9 amplitudes are in good agreement (modulo sign) with the shell-model 
wave functions of Kumar |TJ| . 



The f3is[n] values for 10 B states are given in Table |rVj, and for 10 Be states in Table [V|. The 



neglect of the [321] or lower symmetry states in these nuclei is justified on the grounds that 
this is the leading spatial symmetry of 10 Li, whose isobaric analog state first appears at 21 
MeV in the excitation spectrum of 10 Be. For 10 B, the [42] symmetry states are dominant, and 
addition of the [33] and [411] symmetries improves the energies by only 0.2 MeV. However, 
for 10 Be, these extra components can give significant additional binding of up to 3 MeV. In 
this case, the extra states correspond to low-lying excitations that may not be filtered out 
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of the trial function by a GFMC propagation to limited r, as discussed below. Thus it is 
crucial to carry out the diagonalization in the trial function to get an optimal starting point 
for the GFMC calculation. 

The A = 10 nuclei are exactly midway through the p-shell and are unique in having 
two linearly independent states of the same spatial symmetry contribute, i.e., two X D[42] 
states in 10 Be and two 3 D[42] states in 10 B. To uniquely identify these two possible com- 
binations, we diagonalize Jastrow trial functions containing just the two 25+1 D[42] states 
in the quadrupole moment operator, so the first Pls[u] amplitude reported in Tables [TV|-[V| 
for each state corresponds to the lower (negative) quadrupole eigenvalue and the second to 
the higher (positive) quadrupole eigenvalue. Interestingly, we see that most states where 
these two symmetries contribute are dominated by either one amplitude or the other, e.g., 
the first 3 + in 10 B is almost pure positive quadrupole while the second 3 + is almost pure 
negative quadrupole in composition. Similarly the first 2 + in 10 Be is almost pure negative 
quadrupole while the second is pure positive quadrupole. This is analogous to using the 
LKl labeling scheme applied to 10 B by Kurath in a traditional shell model calculation [19 . 
and there is good agreement as to the dominant [42] symmetry amplitudes with his Table 3. 
Only the first 2 + state in 10 B does not seem to fit in with these observations. 

Energy expectation values with the full \Ev of Eq.(|8|) are evaluated by a Metropolis Monte 
Carlo algorithm as described in . The full wave function at any given spatial configuration 
R = ri, T2, ...ta can be represented by a vector of 2 A x I(A, T) complex numbers, 

v&(R) = 5> a (R)|a> , (18) 

a 

where the ij} a (R) are the complex coefficients of each state |a) which has specific third com- 
ponents of spin and linear combinations of good isospin. For the nuclei considered here, 
this leads to vectors ranging in length from 13,824 for 9 He to 92,160 for 10 Be, although a 
savings of a factor of two is possible in even- /I nuclei by computing in the M = substates. 
The spin, isospin, and tensor operators Of~ 2 ' 6 contained in the potential and other opera- 
tors of interest are sparse matrices in this basis and thus easily evaluated. Kinetic energy 
and spin-orbit operators require the computation of first derivatives and diagonal second 
derivatives of the wave function. These are obtained by evaluating the wave function at 6A 
slightly shifted positions of the coordinates R and taking finite differences. Terms quadratic 
in L require mixed second derivatives and additional wave-function evaluations and finite 
differences. 



B. Green's Function Monte Carlo 

The GFMC method [[H], [21) projects out the exact lowest-energy state, ^o, f° r a given set 
of quantum numbers, using \l/ = lim^oo exp[— (H— E )t]^/ t , where is an optimized trial 
function from the VMC calculation. If the maximum r actually used is large enough, the 
eigenvalue Eq is calculated exactly while other expectation values are generally calculated 
neglecting terms of order — ^t| 2 and higher ||. In contrast, the error in the variational 
energy is of order |^o — ^t\ 2 , and other expectation values calculated with ty? have errors of 
order ~~ ^t|- In the following we present a brief overview of the nuclear GFMC method; 
much more detail may be found in Refs. |2|. [|. 
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We start with a defined as (see Eq.(3.13) of Ref. ||) 

\*T)=SY[(l + U ij +j:U™ I )\* J ), (19) 

i<j k^i,j 

and define the propagated wave function \I/(t) 

tt( r ) = e -(H~E )r^ T = j e -(ff-Bo)Arj» ^ ? (2Q) 

where we have introduced a small time step, Ar = r/n; obviously \l/(r = 0) = and 
\I/(t — > oo) = \P . The \I/(t) is represented by a vector function of R using Eq. fll8|) , and the 
Green's function, G aj g(R, R') is a matrix function of R and R' in spin-isospin space, defined 
as 

G Q/3 (R, R') = (R, a|e-^-^ AT |R', (3) . (21) 

It is calculated with leading errors of order (Ar) 3 . Omitting the spin-isospin indices a, (3 
for brevity, \I/(R n , r) is given by 

tt(R»,r) = J Gf(R n ,R n _ 1 )...G(R 1 ,R )* T (R )dRn-i---dRidRo , (22) 

with the integral being evaluated stochastically. The short-time propagator is constructed 
with the exact two-body propagator and additional terms coming from the three-body in- 
teraction. 

Quantities of interest are evaluated in terms of a "mixed" expectation value between v&y 
and #(t): 



(O) 



Mixed 





OMt)) 




¥(r)> 



_ I rfP»4(Rn)0G(Rn, Rn-l) " " " P(Rl, Rq)^t(Rq) ^ 

/dP n ^(R n )G(R n ,R n _ 1 )---G(R 1 ,R )^ T (Ro) ' 

where P n = R , Ri, • • ■ , R n denotes the 'path', and dP n = c/RodRx • ■ • dR n with the integral 
over the paths being carried out stochastically. The desired expectation values would, of 
course, have ^(r) on both sides; by writing \&(t) = + <5\I/(t) and neglecting terms of 
order [5\l/(r)] 2 , we obtain the approximate expression 

(0(r)> = \^)\i(T)} W <°( r )> Mi - d + KO(r))Mixed - (OH , (24) 

where (O)t is the variational expectation value. 

A special case is the expectation value of the Hamiltonian. The (H(r)) Mixed can be 
re-expressed as 

(y T \e-( H - E ^ T / 2 He- {H - Eo)T/2 \ty T ) 
(#(r)>Mi X ed = ^ T | e -(/f-£oW2 e -(H-£ W2|^ r) ^ E ° > ( 25 ) 

since the propagator exp[— (if — -Eo)t] commutes with the Hamiltonian. Thus (if (T)) Mixcd 
approaches Eq in the limit r — > oo, and the expectation value obeys the variational principle 
for all r. 
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As mentioned above, the propagation is actually made with a simplified Hamiltonian H' 
that includes the CI kinetic energy operator, the AV8' two-nucleon interaction including 
isoscalar Coulomb, and a slightly altered V^ k three-nucleon interaction. The AV8' is a little 
more attractive than AV18, so a slightly more repulsive V^ k is used to keep (H') « (H). 
This ensures the GFMC algorithm will not propagate to excessively large densities due to 
overbinding. Consequently, the upper bound property applies to (H'(r)), and (H — H') 
must be evaluated perturbatively. 

Another complication that arises in the GFMC algorithm is the "fermion sign problem" . 
This arises from the stochastic evaluation of the matrix elements in Eq. fl2"5|) . The Monte 
Carlo techniques used to calculate the path integrals leading to \&(R n ,r) involve only local 
properties, while antisymmetry is a global property. Thus the propagation can mix in 
the boson solution. This has a (much) lower energy than the fermion solution and thus 
is exponentially amplified in subsequent propagations. In the final integration with the 
antisymmetric $?t, the desired fermionic part is projected out, but in the presence of large 
statistical errors that grow exponentially with r. Because the number of pairs that can be 
exchanged grows with A, the sign problem also grows exponentially with increasing A. For 
A>8, the errors grow so fast that convergence in r cannot be achieved. 

To remedy this situation, a "constrained path" approximation has been developed 0. 
The basic idea of the constrained-path method is to discard those configurations that, in 
future generations, will contribute only noise to expectation values. Many tests of the con- 
strained path have been made and it usually gives results that are consistent with uncon- 
strained propagation, within statistical errors, although there are cases in which it converges 
to the wrong energy 0]. Up to now, this problem was always solved by using a small number, 
n u = 10 to 20, of unconstrained steps before evaluating expectation values. These few un- 
constrained steps, out of 400 or more total steps, were enough to damp out errors introduced 
by the constraint in all the tests that we had done. The statistical errors in calculations 
with n u = 10 are not substantially greater than for n u = 0. However for A = 10, the error 
from n u = 20 is 60-90% larger than that from n u = 10, i.e., 2.5-3.5 more configurations 
must be used to get the same error. In the present work, we find that A = 9, 10 calcula- 
tions with the AV18/IL3 Hamiltonian (which has the strongest long-range part of Vijk) are 
particularly sensitive to n u , and most of these are made with n u = 20. Tests made with the 
other Hamiltonians show that n u = 10 is enough, except for 9 ' 10 He. 

Previously we had found that obtaining reliable results for 8 He was particularly dif- 
ficult. In the present work we find that this persists for 9 ' 10 He which show large statistical 
fluctuations. In addition the constrained-path method appears to be less reliable than for 
other nuclei. For example, we made a calculation for 10 He with AV18/IL2 and a that 
contained no NNN correlation. Even though n u = 20 was used, the result was overbound 
by 1.5(6) MeV. (In Ref. || a calculation of 6 Li with a much worse *&t, which contained 
no NN tensor correlations, also gave an overbound result but this was corrected with only 
n u = 10.) The statistical errors using n u = 40 are very large, making a confirmation of the 
n u = 20 results presented here difficult; we have attempted this for AV18/IL2 and obtain 
no significant change in the binding energy within a statistical error of 0.7 MeV. 

Figure |^ shows the progress with increasing r of typical constrained GFMC calculations, 
in this case for various states of 10 B using the AV18/IL2 Hamiltonian. The values shown at 
t = are the VMC values using ^t- The GFMC very rapidly makes a large improvement 
on these energies; by r = 0.01 MeV -1 , the ^>t energies have been reduced by ~20 MeV. This 
rapid improvement corresponds to the removal of small admixtures of states with excitation 
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energies ~ 1 GeV from ^t- Typically, averages over the r > 0.1 MeV -1 values are used as 
the GFMC energy. The standard deviation, computed using block averaging, of all of the 
individual energies for these r values is used to compute the corresponding statistical error. 
The solid lines show these averages; the corresponding dashed lines show the statistical 
errors. 



V. ENERGY RESULTS 

In this section we present GFMC energy results in A = 9, 10 nuclei for the simplified AV8' 
NN Hamiltonian, the full AV18 interaction, and AV18 with each of the NNN interactions 
UIX, IL2, IL3, and IL4. (Results for A < 8 nuclei may be found in Ref. 0].) The total 
energies for 18 states are shown in Table |VT] and in Fig. |]; excitation spectra are shown 
in Fig. |5]. Not all states have been calculated with all the Hamiltonians. Numbers in 
parentheses are the Monte Carlo statistical errors; in addition there may be systematic 
errors from the GFMC algorithm of the order of 1-2% percent as discussed above. 

The energies of the particle-stable ground states calculated with just the NN interactions 
are underbound by 8-12 MeV in A = 9 and 12-16 MeV in A = 10, with the AV8' being 
somewhat more attractive than AV18. Addition of the older UIX NNN interaction to AV18 
reduces this discrepancy significantly to 2-4 MeV in A = 9 and 4-6 MeV in A = 10. 
However, addition of the IL2-IL4 models to AV18 results in much better energies, some 
high and some low, but generally within 2 MeV. In particular, the AV18/IL2 and AV18/IL3 
Hamiltonians come within ±3% of the experimental bindings for the ground states of 9 Li, 
9 Be, 10 Be, and 10 B, which is about as good as can be expected, given the 1-2% total errors 
in the GFMC calculation. 

Despite their systematic underbinding, the older Hamiltonians were able to give the 
correct ordering of excited states in the A < 8 nuclei, albeit usually with insufficient splitting 
between states that are spin-orbit partners ||. In the A = 9, 10 nuclei this no longer holds 
true. The most dramatic case is the question of the proper ground state for 10 B: AV8', AV18, 
and AV18/UIX all clearly and wrongly predict the (1 + ; 0) state to be lower than the (3 + ; 0) 
state, while all the AV18/IL2-IL4 models correctly reproduce the experimental ordering. 
This result for the AV18 interaction was first obtained by Navratil who also finds it for the 



CD-Bonn interaction ||22| . This incorrect ordering of states by the two-body Hamiltonians 
also seems to be present in the case of the (| ; |) and (| ; |) excitations of 9 Be, and may 
be a problem for the proper ground state of 9 Li, although with the statistical errors of the 
calculation it is not so clear. The AV18/UIX and AV18/IL2-IL4 models obtain these states 
in the proper order, with one exception, and generally with fairly good spin-orbit splittings. 

As noted in Sec. II, two narrow, particle-unstable states have been observed in 9 Li 
without definite spin assignments. On the basis of the close agreement between the AV18/IL2 
energies for the lowest (| ; |) and (| ; \) states and the observed excitations, as illustrated 
in Fig. |5], we think these are in fact the most likely assignments for these states. Also, our 
one calculation of the expected (3 + ; 1) state in 10 Be lines up with the observed state at 9.4 
MeV, suggesting that is its proper spin assignment. 

The current calculations significantly underbind 9,10 He, even when the AV18/IL2-IL4 
models are used. Because both these nuclei are particle unstable, they should properly be 
calculated as resonant states, so the current GFMC results may be somewhat less reliable 
than for the other ground states. However, the experimental widths of these states are not 
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large, and our calculations for other comparably narrow states have generally been good. 
It may be that some significant admixture of lsj/2 orbitals from the sd-shell should be 
introduced in the Jastrow trial function, as suggested by recent shell- model studies [^3|; we 
have not yet attempted such a modification. In addition, as was discussed in the previous 
section, the constrained-path propagation may not be as reliable for these nuclei. If the 
present results for 9,10 He are correct, then there is a need for further tuning of the isospin 
dependence of the NNN interactions. 

In addition to the total energies discussed above, we have calculated the energy differences 
between isobaric analog states perturbatively by evaluating the expectation values of the 
electromagnetic and strong CSB and CD parts of the AV18 Hamiltonian in the wave function 
of the T z = — T nucleus. The CSB and CD terms can induce corresponding changes in the 
nuclear wave functions, leading to higher-order perturbative corrections to the splitting of 
isospin mass multiplets, but it is difficult for us to estimate these higher-order effects reliably 
in either VMC or GFMC calculations. The energies for a given isomultiplet of states can be 
expanded as 

E a ,t(T z ) = E a n (A,T)Q n (T,T z ) , (26) 

n<2T 

where Q = 1, Q 1 = T z , and Q 2 = |(3T 2 2 - T 2 ) are isoscalar, isovector, and isotensor 
terms The isovector coefficients ai(A,T) are shown in Table |V1I| , while the isotensor 
coefficients az(A,T), are given in Table [VIII| . 

The dominant contribution in each case is from the regular Coulomb interaction between 
protons, but there are additional contributions of ~ 5% coming from the other electromag- 
netic and strong terms. Because the AV18 and AV18/UIX models underbind these nuclei, 
the density distributions tend to be more diffuse and the Coulomb interaction is reduced, 
leaving a sizable discrepancy with experiment for the isovector coefficients. The AV18/IL2- 
IL4 models give much better binding energies and consequently larger isovector coefficients 
which are within 5% of the experimental values. However, the isotensor coefficients come 
out a bit too large. 



VI. MOMENTS AND DENSITY DISTRIBUTIONS 



Point proton and neutron rms radii are shown in Tables [IX| and |X| The AV18/IL2-IL4 



Hamiltonians give significantly smaller radii than the older models, presumably because of 
their greater binding. The experimentally measured charge radii PB| for the 9 Be and 10 B 
ground states are in good agreement with the calculations from the new models. 

Quadrupole moments are shown in Table DHL as evaluated in impulse approximation. 



Experimental values are from Refs. || [7J. The contribution of two-body charge operators to 
the quadrupole moment are expected to be small, as in the deuteron and 6 Li [f|, |2"7| , because 
both the isoscalar and isovector charge operators are relativistic corrections of order (v/c) 2 . 
Expectation values of the quadrupole operator tend to have a much larger variance than 
radii because of the cancellations from the Y 2 o operator. Given this caveat, the agreement 
with the measured quadrupole moments is reasonable. 

The impulse approximation magnetic moments are shown in Table |XII| ; experimental val- 



ues are from Refs. || [7| |31|. Experience in lighter nuclei shows that there are significant 
two-body current contributions to the magnetic operators, particularly the isovector part. 
In the trinucleons, the isoscalar portion is boosted by +0.034/^ or « 8%, while the isovector 



portion is corrected by — 0.778/iAr or « 18% p8| . However, in the isoscalar 6 Li case, the 
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correction is a tiny +0.003 /ijy [27 ■ On this basis, we expect substantial improvement in the 
9 Li and 9 C results shown in Table Xll when two-body contributions to the current operator 
are added, while the already good 10 B results will not be significantly affected. However, the 
reasonable agreement between the IA calculations and experiment for 9 Be may not survive. 

The nucleon densities for 9 Li and 9 Be are shown in Fig. on both a linear (left) and 
logarithmic (right) scale. The densities are normalized such that the integrated value equals 
the appropriate total value of N or Z. These densities are rather similar to those calculated 
earlier for 8 Li and 8 Be with the former more peaked at the origin, and the latter showing 
a broad flat interior. The protons in 9 Li cluster at the origin because two out of three are 
confined to the a core, while the neutrons show a substantial skin with a hint of a p-shell 
peak around 1 fm. The low central densities for both protons and neutrons in 9 Be are 
probably due to the large 2a component in its 8 Be core. The densities for 10 Be and 10 B 
are shown in Fig. |7|. They have somewhat higher central values than in 9 Be, but are also 
broader and flatter than 9 Li, presumably due to the dominant 8 Be core. 



VII. CONCLUSIONS 

We have made quantum Monte Carlo calculations for the ground states and low-lying 
excited states of A = 9, 10 nuclei, based on realistic two- and three-nucleon potentials. The 
new Illinois NNN models do a good job of explaining the binding energies of the ten best- 
known A = 9, 10 states of different quantum numbers, but the results for 9 ' 10 He may indicate 
a deficiency in the models. If we also include the 17 narrow states to which these forces were 
fit in A < 8 nuclei [||, we find rms deviations from the experimental energies of 0.6 MeV 
for AV18/IL2, 0.8 MeV for AV18/IL3, and 1.0 MeV for AV18/IL4. This contrasts with rms 
deviations of 3.2 MeV for AV18/UIX and 9.9 MeV for AV18 alone. Thus AV18/IL2 becomes 
our preferred model Hamiltonian with AV18/IL3 almost as good, although AV18/IL4 will 
be somewhat deprecated. 

We conclude that a fairly consistent picture of nuclear binding can be constructed for 
A < 10 nuclei using a single Hamiltonian and a single computational framework. This 
applies also to the energy differences among isobaric multiplets, which are well-reproduced. 
Electromagnetic moments, within the limitation of the impulse approximation, are in fairly 
good agreement with experimental data. 

Significant challenges for the future will be the computation of second or higher excited 
states of the same quantum number, and intruder states of unnatural parity; both of these 
types of states first become particle-stable in the A = 10 nuclei. We also need to evaluate 
electroweak transition rates in these nuclei, as has been done previously in the A = 6, 7 
systems p7| , |2Ti| . The present calculations are at the limit of our current computer resources, 
both in CPU time and resident memory. However, the present codes are in principle capable 
of computing nuclei like n Be and 12 C; such calculations may become feasible with the next 
generation of computers. 
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TABLE I: Permutation symmetry terms for LS'-coupling in A = 9, 10 nuclei and corresponding 
spin states. 

A [n] L (T, S) highest symmetry states 

9 [41] 1,2,3,4 (1,1) 9 Be( i-_ |- } 

[32] 1,2,3 (§,!)(!, §)(^) 9 Li(r-D 
[311] 0,2 (|,|)(|,|)(|,|)(|,|) i 

[221] 1 (5, 2X2' 2X2 ' 2X2' 2X2' 2X2' 2) He(g '2 ) 

10 [42] 0, 2 2 , 3, 4 (1, 0)(0, 1) 10 Be(0+, 2+- 4+), 10 B(1+- 5+) 

[33] 1,3 (1, 1X0,0) 10 Bc(l+) 
[411] 1,3 (1,1)(0,0) 

[321] 1,2 (2,1)(2,0)(1,2)(1,1) 2 (1,0)(0,2)(0,1) 10 Li(0+- 3+), 10 B(0+) 

[222] (3,0)(2,1)(1,2)(1,0)(0,3)(0,1) 10 He(0+) 



TABLE II: VMC diagonalization for Plsm components in Be. 



TIT 
J 


2 P[41] 


2 D[41] 


2 F[41] 


2 G[41] 


4 P[32] 


4 D[32] 


4 F[32] 


2 P[32] 


2 D[32] 2 F[32] 


3- 

2 


0.936 


-0.337 






0.035 


-0.024 


0.047 


-0.026 


0.049 


5~ 
2 




0.952 


0.273 




-0.064 


-0.040 


0.060 




0.062 -0.041 


1- 
2 


0.990 








0.117 


-0.039 




0.002 




7- 
2 






0.868 


0.488 




0.019 


-0.064 




0.034 


3-* 
2 


0.356 


0.921 






0.087 


-0.098 


0.052 -0.050 -0.003 


5 5(S 

2 




-0.248 


0.925 




0.273 


0.007 


0.045 




0.037 -0.051 


9- 
2 








0.978 






0.207 






1-* 
2 


-0.104 








0.905 


0.286 




-0.153 




7-* 
2 






-0.474 


0.836 




-0.260 


0.009 




0.042 




4 S[311] 4 D[311] 2 S[311] 2 D[311] 6 P[221] 4 P[221] 2 P[221] 






3- 
2 


0.046 


-0.025 




0.011 


0.029 


-0.010 


0.010 






5- 
2 




0.053 




-0.005 


-0.025 


-0.015 








1- 
2 




-0.054 


-0.023 






-0.038 


0.024 






7- 
2 




0.057 






-0.015 










3-* 
2 


0.046 


0.035 




-0.000 


-0.005 


0.000 


0.010 






5 5(S 

2 

9- 




0.003 




-0.035 


-0.017 


0.021 








2 

1-* 
2 




-0.112 


0.216 






-0.077 


0.005 






7-* 
2 




-0.080 






-0.015 
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TABLE III: VMC diagonalization for PlsWi] components in 9 Li. 
r 2 P[32] 2 D[32] 2 F[32] 4 S[311] 4 D[311] 2 S[311] 2 D[311] 4 P[221] 2 P[221] 
~Y 0.859 -0.089 -0.446 -0.070 -0.198 0.107 -0.004 

i~ 0.841 -0.487 0.219 0.028 0.085 

§~ 0.968 -0.234 -0.045 -0.060 -0.042 

§~* -0.130 0.844 -0.489 0.016 0.160 -0.017 -0.075 

\~ 0.826 -0.563 

\~* -0.029 0.333 0.934 -0.109 0.054 

f ~* -0.088 -0.518 0.821 0.163 -0.154 



TABLE IV: VMC diagonalization for 0LS[n] components in B. 





3 S[42] 3 D[42] 


3 D[42] 3 F[42] 


3 G[42] 


X P[33] 


^[33] 


X P[411] 


^[411] 


3+ 


0.036 


0.995 -0.086 


0.001 




0.016 




0.003 


1+ 


0.889 -0.225 


0.368 




0.068 




-0.137 




l+* 


0.179 0.922 


0.024 




-0.021 




-0.342 




2+ 


0.467 


0.884 0.023 












3+* 


0.917 


0.368 0.016 


-0.047 




0.021 




0.142 


-^+** 


-0.042 -0.014 


0.998 




-0.025 




-0.038 




2+* 


-0.141 


0.978 -0.154 












4+ 




0.901 


0.433 











TABLE V: VMC diagonalization for /?£,s[„] components in 10 Be. 
,r ^S[42j ^[42] X D[42] X F[42] ^G[42] 3 P[33] 3 F[33] 3 P[411] 3 F[411] 
"5+ 0.812 0.109 0.573 

2+ 0.944 -0.071 -0.050 0.005 0.288 -0.138 

2+* -0.035 0.998 -0.030 0.019 0.032 0.017 

0+* -0.462 -0.492 0.738 

3+ 0.969 0.102 -0.225 

4+ 0.867 -0.228 -0.443 
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TABLE VI: Total GFMC energies in MeV. 





AV8' 


AV18 


UIX 


IL2 


IL3 


IL4 


Expt 


9 He(H 


-22.5(2) 


-20.7(3) 




-28.7(3) 


-28.2(5) 


-28.8(4) 


-30.21(8) 


9 Li(F) 


-36.6(2) 


-33.7(3) 


-40.9(3) 


-46.0(4) 


-46.7(5) 


-47.6(4) 


-45.34 




-36.8(2) 


-34.0(3) 


-39.4(3) 


-43.8(4) 


-44.2(4) 


-44.5(5) 


-42.65 


9 Li(D 


-35.0(2) 


-32.1(3) 


-37.9(4) 


-41.1(4) 


-41.0(4) 


-40.6(4) 


-39.96(?) 


9 Li(D 


-32.0(2) 


-29.7(3) 


-35.2(3) 


-39.0(4) 


-38.7(4) 


-40.8(4) 


-38.91(?) 


9 Be(r) 


-49.9(2) 


-46.4(4) 


-55.1(3) 


-58.2(5) 


-57.8(5) 


-58.0(6) 


-58.16 


9 Be(£~) 


-47.8(2) 


-43.5(3) 


-51.3(5) 


-55.8(5) 


-54.1(4) 


-55.1(5) 


-55.73 


ot ^2 > 


—48.2(2) 


—45.0(4) 


—50 9(6) 


-54.3(4) 


—55.7(5) 


—55.5(4) 


—55.36 


ae \2 > 


— ^tO.Oy^ ) 










— O-L.Ol / ) 


— oi . to 


Be(2 J 
rie^U J 


/in i f o\ 
— 4U.1^Z J 


— OD. l{o) 




A1 Cii R\ 

— 4 1 .\)yo ) 




—4 / .0{0 ) 




— ZL.Z[Z ) 


in 

— iy.o(o j 




— Zo.4(o J 


— Z 1 .Z[0 ) 


07 7(K\ 
— Z t .i[0) 


on qa 
— oU.o4 


oe^U ) 


— DO. L[Z ) 


— OZ.U(0 J 


— o\).lyo ) 


— DD.8( 1 ) 


— 00. Z( 1 ) 


A { R\ 
— O / .4^0 J 


dA OQ 

— D4.yo 


10 Be(2+) 


-51.9(2) 


-47.7(5) 


-57.1(6) 


-61.8(5) 


-59.9(6) 


-61.1(6) 


-61.61 


10 Be(3+) 








-55.7(6) 






-55.58 


10 B(3+) 


-53.2(3) 


-48.6(6) 


-59.0(4) 


-65.6(5) 


-64.1(5) 


-65.6(6) 


-64.75 


10 B(1+) 


-55.7(3) 


-51.6(6) 


-60.3(5) 


-64.7(4) 


-62.8(5) 


-65.1(5) 


-64.03 


10 B(2+) 


-52.2(2) 


-47.2(5) 




-61.7(5) 


-60.6(7) 


-62.4(5) 


-61.16 


io B(4 +) 


-50.0(3) 


-45.0(5) 




-60.0(5) 


-58.2(6) 


-61.5(4) 


-58.72 



TABLE VII: Isovector energies in keV. 

AV18 UIX IL2 IL3 IL4 EXPT 

9 Li(§~) 1900(8) 1860(7) 2000(7) 2041(7) 2113(7) 2104 

9 Li(l~) 1831(8) 1804(8) 2003(7) 1891(7) 2047(8) 1946 

9 Be(f~) 1751(10) 1738(9) 1857(10) 1833(10) 1878(14) 1851 

9 Be(|") 1700(10) 1704(8) 1753(12) 1754(12) 1869(13) 1783 

10 Be(0 + ) 1970(10) 2115(12) 2235(12) 2140(13) 2307(11) 2329 

10 Be(2+) 2111(11) 1970(12) 2189(12) 2124(11) 2079(11) 2322 



TABLE VIII: Isotensor energies in keV. 
AV18 UIX IL2 IL3 IL4 EXPT 

9 Li(f~) 182(9) 182(9) 197(8) 202(11) 204(9) 176 
9 Li(i~) 172(7) 181(8) 193(8) 181(8) 184(10) 160 
10 Be(0+) 257(25) 205(33) 250(24) 284(29) 314(27) 241 
10 Be(2+) 223(27) 196(35) 216(18) 183(23) 208(28) 199 
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TABLE IX: Proton rms radii in fm. 





AV8' 


UIX 


IL2 


IL3 


IL4 


EXPT 


9 Li(F) 


2.19(1) 


2.20(1) 


2.04(1) 


2.02(1) 


1.94(1) 




9 Li(i~) 


2.27(1) 


2.23(1) 


2.11(1) 


2.29(1) 


2.07(1) 




9 Be(r) 


2.41(1) 


2.41(1) 


2.38(1) 


2.36(1) 


2.33(1) 


2.40(1) 


uc l 2 ) 


2.41(1) 


2.41(1) 


2.38(1) 


2.46(1) 


2.27(1) 




10 Be(0+) 


2.38(1) 


2.30(1) 


2.21(1) 


2.32(1) 


2.19(1) 




10 Be(2+) 


2.33(1) 


2.42(1) 


2.26(1) 


2.29(1) 


2.28(1) 




10 B(3+) 


2.40(1) 


2.49(1) 


2.33(1) 


2.38(1) 


2.27(1) 


2.33(12) 


10 B(1+) 


2.45(1) 


2.48(1) 


2.49(1) 


2.52(1) 


2.42(1) 




10 B(2+) 


2.46(1) 




2.38(1) 


2.59(1) 


2.29(1) 




io B(4 +) 


2.42(1) 




2.30(0) 


2.41(1) 


2.15(0) 





TABLE X: Neutron rms radii in fm. 

AV18 UIX IL2 IL3 IL4 



9 Li(F) 

9 Li(D 
9 Be(|") 

9 Be(|") 

10 Be(0+) 

10 Be(2+) 

10 B(3+) 

10 B(1+) 

10 B(2+) 
io B ( 4 +) 



2.72(1) 
2.84(1) 
2.63(1) 
2.67(1) 
2.76(1) 
2.66(1) 
2.40(1) 
2.45(1) 
2.46(1) 
2.42(1) 



2.76(1) 
2.80(1) 
2.63(1) 
2.65(1) 
2.56(1) 
2.78(1) 
2.49(1) 
2.48(1) 



2.57(1) 
2.56(1) 
2.55(1) 
2.62(1) 
2.47(1) 
2.55(1) 
2.33(1) 
2.49(1) 
2.38(1) 
2.30(0) 



2.52(1 
2.76(1 
2.56(1 
2.65(1 
2.59(1 
2.60(1 
2.38(1 
2.52(1 
2.59(1 
2.41(1 



2.39(1 
2.53(1 
2.52(1 
2.48(1 
2.45(1 
2.62(1 
2.27(1 
2.42(1 
2.29(1 
2.15(0 





TABLE 


XI: Quadrupole moments in fm 2 . 






AV8' 


UIX 


IL2 


IL3 


IL4 


EXPT 


9 Li(§-) 


-3.1(1) 


-2.9(1) 


-2.7(1) 


-2.7(1) 


-2.5(1) 


-2.7(1) 


9 Be(|-) 


5.0(3) 


7.4(2) 


8.5(3) 


5.7(2) 


6.8(3) 


5.9(1) 


9 Be(D 


-2.3(2) 


-2.6(2) 


-2.9(2) 


-4.8(2) 


-4.1(2) 




9 C(1") 


-3.7(3) 


-5.5(3) 


-5.1(3) 


-4.4(3) 


-5.3(3) 




10 Be(2+) 


-10.1(3) 


-8.0(4) 


-5.0(4) 


-9.6(4) 


7.2(2) 




10 B(3+) 


8.7(2) 


12.0(2) 


9.5(2) 


9.4(2) 


8.8(2) 


8.5(1) 


10 B(1+) 


1.2(1) 


3.2(2) 


3.1(1) 


3.2(1) 


2.1(1) 




10 B(2+) 


-1.9(2) 




-0.6(2) 


-1.8(3) 


-1.8(2) 




io B ( 4 +) 


3.9(3) 




3.8(2) 


5.4(4) 


3.8(2) 
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TABLE XII: Magnetic moments in /ijy. 





AV8' 


UIX 


IL2 


IL3 


IL4 


EXPT 


9 LKr) 


2.91(1) 


2.54(2) 


2.47(2) 


2.54(2) 


2.54(2) 


3.44(0) 


9 Li(I-) 


-0.23(2) 


-0.19(3) 


-0.29(2) 


-0.25(3) 


-0.24(3) 




"Bed") 


-1.35(2) 


-1.11(1) 


-1.14(2) 


-1.03(2) 


-1.17(3) 


-1.18(0) 


9 Be(f-) 


-0.95(1) 


-0.96(1) 


-0.88(1) 


-0.84(1) 


-0.85(1) 




9 c(|") 


-1.08(3) 


-0.71(5) 


-0.63(4) 


-0.72(4) 


-0.70(4) 


-1.39(0) 


10 B(3+) 


1.82(1) 




1.80(1) 




1.80(1) 


1.80(0) 


10 B(1+) 


0.75(1) 


0.78(1) 


0.74(1) 




0.77(1) 


0.63(12) 
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Pieper, Varga, & Wiringa — Fig. 1 
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FIG. 1: First natural-parity states in the experimental spectra of A = 9 nuclei. 
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Pieper, Varga, & Wiringa — Fig. 2 
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FIG. 2: Low-lying natural-parity states in the experimental spectra of A = 10 nuclei. 
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Pieper, Varga & Wiringa ~ Fig. 3 
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FIG. 3: GFMC E(t) for states of 10 B calculated for the AV18/IL2 Hamiltonian. The solid lines 
show the averages that are reported in the text; the dashed lines show the corresponding statistical 
errors. 
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Pieper, Varga & Wiringa - Fig. 6 
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FIG. 6: Proton and neutron densities for A = 9 nuclei on both linear (left) and logarithmic ( 
scales. 
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Pieper, Varga & Wiringa - Fig. 7 
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FIG. 7: Proton and neutron densities for A = 10 nuclei on both linear (left) and logarithmic (: 
scales. 
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